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ABSTRACT 

The paper describes the application and assessment of the 
recently-developed CAP-TSD transonic small-disturbance code 
for flutter prediction. The CAP-TSD code has been developed for 
aeroelastic analysis of complete aircraft configurations and was 
previously applied to the calculation of steady and unsteady 
pressures with favorable results. Generalized aerodynamic 
forces and flutter characteristics are calculated and compared in 
the present study with linear theory results and with 
experimental data for a 45° sweptback wing. These results are 
in good agreement with the experimental flutter data which is 
the first step toward validating CAP-TSD for general transonic 
aeroelastic applications. The paper presents these results and 
comparisons along with general remarks regarding modern wing 
flutter analysis by computational fluid dynamics methods. 


NOMENCLATURE 


Aij 

generalized aerodynamic force coefficient from 
the surface integral of hi Apj/lc 2 pU 2 /2J 

bo 

reference length, c/2 

c 

reference length; root chord 


coefficient of lifting pressure, Ap/(pU 2 /2) 

k 

reduced frequency, co bo/U 

M 

freestream Mach number 

Ap, Apj 

lifting pressure; time marching value and that 
for mode hj, respectively; positive up 

qi(t) 

generalized coordinate of motion for mode i 

t 

time 

u 

freestream speed 

x, y. z 

right-hand orthogonal coordinates 

P 

ratio of wing mass to mass of air in the truncated 
cone that encloses the wing 

P 

freestream flow density 

(0, coj, (Dot 

circular frequency, in general; of mode j and of 
the first torsion mode, respectively 


INTRODUCTION 

Research on the application of computational fluid dynamics 
(CFD) methods to unsteady flows and aeroelastic analysis is 
presently being actively pursued. Edwards and Thomas (1) gave 
a recent survey, for example, on computational methods for 
unsteady transonic flows with emphasis on applications to 
aeroelastic analysis and flutter prediction. The transonic speed 
range is of primary interest because the flutter dynamic 
pressure is typically critical (i.e., lower) there. The main 
effort, especially for three-dimensional configurations, has 
been at the transonic small disturbance (TSD) equation level, of 
which the XTRAN3S program is an important example (2). For 
the higher equation levels such as the Euler and Navier-Stokes 
equations, efforts on aeroelastic applications have been limited 
to simple two-dimensional airfoils because of the larger 
computational cost involved. Two recent efforts are reported by 
Bendiksen and Kousen (3) and by Wu, Kaza, and Sankar (4). 

The advantage of the TSD formulation, though, is the 
relatively low computational cost, the simplicity of the gridding 
and geometry preprocessing, and the ability to treat complete 
aircraft configurations. The XTRAN3S code, for example, 
employs an alternating-direction implicit (ADI) finite- 
difference algorithm for solution of the TSD equation, with 
several terms treated explicitly. This type of solution has a 
numerical stability restriction, however, that requires a large 
number of sufficiently small time steps, often much smaller 
than are needed for a time-accurate solution of the physical 
process. Batina (5) described the development of a 
time-accurate approximate factorization (AF) algorithm applied 
to the TSD equation. The AF algorithm was developed to alleviate 
or avoid the numerical stability restriction of the ADI 
algorithm. Furthermore, the rate of convergence per time step 
is greatly enhanced by the AF procedure, which reduces 
proportionately the cost of computation. 

The AF algorithm has subsequently been developed Into a new 
computer code called CAP-TSD (for Computational Aeroelasticlty 
Erogram - Lransonic Email {disturbance) for transonic 
aeroelastic analysis of complete aircraft configurations (6). 



The CAP-TSO program emphasizes a combination of economy, 
stability, and accuracy of calculation. CAP-TSD has been used to 
calculate steady and unsteady pressures on wings and 
configurations at subsonic, transonic, and supersonic Mach 
numbers. Comparisons of these results with other methods and 
with experimental data have been favorable (Refs. (6) and 
(7)). However, the CAP-TSD code has been developed primarily 
for aeroelastic analysis. Such analysis involves the coupling of 
the aerodynamics with the structural characteristics of the 
configuration under consideration. The resulting equations of 
motion for a time-domain or time-marching aeroelastic 
analysis are based upon the aircraft natural vibration modes. 
These equations are integrated in time along with the finite- 
difference solution of the flow field. Initial conditions for each 
mode are input and free decay transients are calculated. 
Aeroelastic stability is then deduced from the free decay records 
or time histories. Both the underlying theory and the numerical 
procedures require evaluation. Thus, the purpose of the present 
paper is to report on the preliminary results of this evaluation. 
Generalized aerodynamic forces and flutter boundaries are 
presented for a 45° sweptback wing. Comparisons of these 
results with parallel linear theory calculations as well as with 
the experimental flutter data of Yates, et al. (8) provide an 
assessment of CAP-TSD for aeroelastic applications. The paper 
presents these results and comparisons along with general 
remarks regarding modern wing flutter analysis by CFD 
methods. 


COMPUTATIONAL PROCEDURES 

In this section, the computational procedures are described 
including the CAP-TSD code, the aeroelastic equations of motion, 
the time-marching solution of these equations, and the modal 
identification of the resulting free decay transients. Although 
the emphasis of the discussion is on the CAP-TSD code, the 
general computational procedures regarding time-marching 
flutter analysis apply to other CFD codes as well. 


C Aff-T SP. C adfl 

The CAP-TSD code is a finite-difference program which 
solves the general-frequency modified transonic small- 
disturbance (TSD) equation. The TSD potential equation is 
defined by 

M 2 (♦, + 2<> x ) t = ((1 - M\ + F^+ G$ x 


+ (* y + H4Mt> y ) y + (4> r ) z (i) 


Several choices are available for the coefficients F, G, and H 
depending upon the assumptions used in deriving the TSD 
equation. For transonic applications, the coefficients are herein 
defined as 

F = -^(y+l)M 2 (2a) 


G = i(y-3)M 2 (2b) 


H = -(y- 1)M J (2c) 


The linear potential equation is recovered by simply setting F, 
G, and H equal to zero. 

Equation (1) is solved within CAP-TSD by a time-accurate 
approximate factorization (AF) algorithm developed by Batina 
(5). In Refs. (5) to (7), the AF algorithm was shown to be 
efficient for application to steady or unsteady transonic flow 
problems. It can provide accurate solutions in only several 
hundred time steps yielding a significant computational cost 
savings when compared to alternative methods. Recently several 
algorithm modifications have been made which improve the 
stability of the AF algorithm and the accuracy of the results (9). 
These algorithm modifications include: (1) Engquist-Osher 

(E-O) type-dependent differencing to more accurately and 
efficiently treat regions of supersonic flow, (2) extension of the 
E-0 switch for second-order-accurate upwind differencing in 
supersonic regions to improve the accuracy of the results, (3) 
nonreflecting far field boundary conditions for more accurate 
unsteady applications, and (4) several modifications which 
accelerate convergence to steady state. The capabilities 
employed in the present study include the E-0 switch and the 
nonreflecting boundary conditions. The CAP-TSD code can treat 
configurations with arbitrary combinations of lifting surfaces 
and bodies including canard, wing, tail, control surfaces, tip 
launchers, pylons, fuselage, stores, and nacelles. Results have 
been presented for several complex aircraft configurations in 
Ref. (6). The calculated results were in good agreement with 
available experimental pressure data which validated CAP-TSD 
for multiple component applications with mutual aerodynamic 
interference effects. 


Equations of Motion 

The aeroelastic equations of motion which have been 
incorporated in CAP-TSD are based on a right-hand orthogonal 
coordinate system with the x-direction defined as positive 
downstream and the z-direction positive upward. The 
presentation herein is limited to the case of an Isolated wing 
with motion in the z-direction from an undisturbed position in 
the z - 0 plane. The general motion of the wing is assumed to be 
described by the separation of time and space variables 


z(x, y, t) « f(t) Z(x,y) ( 3 ) 


and further that the motion can be well approximated as a finite 
modal series 

z(x, y, t) - ^ q. (t) hj (x, y) ( 4 ) 


where for each mode, hi is the mode shape and qi is the 
generalized coordinate of modal motion. The equations of motion 
are then formulated by considering Lagrange's equations. The 
principal of virtual work as expressed by a set of Lagrange's 
equations (one equation per mode i) is 


2 



(5) 


( 12 ) 


d / 3T 3T dU 


Qi 


a -l. 2.-0 


Qi 


PU 2 2f. Ap dS 
2 j ' pU 2 /2 c 2 


where T is the kinetic energy, U is the potential energy of the 
structure, and Qi is the generalized force associated with qi, 
which is obtained from the virtual work due to the external 
(aerodynamic) forces. The kinetic energy is 

(.) 

s * 


where S is the planform area and dm/dS is the area distribution 
of wing mass. The potential energy is 

u “ J i ° (x,y) h i < x *y » 2 dS < 7 ) 

s ' 


is the generalized aerodynamic force associated with the lifting 
pressure Ap, and with the virtual woik weighting mode h|. The 
resulting set of equations of motion is 

Mq + Kq = Q 

(13) 

q T = [q, q 2 •••] 

where M is the generalized mass matrix, K is the stiffness 
matrix, and Q is the vector of generalized aerodynamic forces. 
From Eq. (13) 

q * - M’ 1 K q + M' 1 Q (14) 


where o(x,y) represents the effective stiffness rate of the wing 
elastic restoring force. 

Application of Lagrange's equations in the absence of external 
forces gives 

£ i vJ^"i dSt 2 < iiH h J dS * 0 m 

is is 


Time-Marchinn Aeroelastic Solution 

The aeroelastic solution procedure implemented within CAP- 
TSD for integrating Eq. (14) is similar to that described by 
Edwards, et. at. (10). Reference (10) describes for a two- 
dimensional, two-degree-of-freedom system an aeroelastic 
solution in terms of a state equation formulation. Here, by a 
parallel formulation, a linear state equation is developed from 
Eq. (14). Each element of Eq. (14) is a normal mode equation 
which may be expressed in first-order state-space form as 


For a linear (or nearly linear) structure in the absence of 
external forces 


Xj = Ax. + Buj 


(15) 


- MM 

■qje . 



(9) 


where to is the circular frequency of vibratory motion. From 
Eq. (8) this gives for each normal mode j 

" 0 > 

is is 

which is the usual Rayleigh-type representation for replacing 
4 the elastic restoring forces involving o in terms of the mass- 

inertia terms involving mp dm/dS. 

t 

The virtual work 8W done by the lifting pressure Ap acting 
through the virtual displacement hi 8 qi is then 


where 

*i * tqi Qil T 


and 


r° 


B . m -tpU 4 2 
B-m. — c 



u j “ [ J AC p hj dS/c 2 ] 
s 


5W = jAph.8qjdS«^ (qj Q,) 5 qj (ii) 

s * 


AC. 


Ap 

pU 2 /2 


from which 
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In these definitions, mi and ki are elements of the mass and 
stiffness matrices, respectively, corresponding to mode i. 
Equation (15) is a finite-dimensional linear differential 
equation and its solution is given by 


t 

x.(t) = O(t) x.(0) + J exp [ A (t - T)] Bu (T) dl (16) 

O 


The state transition matrix <t>(t) - exp [At], in general, can be 
calculated to any assigned accuracy by using a sufficient number 
of terms of the series expansion of the matrix exponential 
function. For the aeroelastic problem considered here, <t*(t) is 
computed exactly using simple closed-form expressions for each 
element of the matrix. As explained in Ref. (10), the first term 
in Eq. (16) is the homogeneous response portion of Eq. (15), 
while the second term is a convolution integral which represents 
the forced response. Numerically, the solution is advanced from 
any time step n to step n + 1 , by 

x. [(n + l)At] = C> (At) x. (nAt) 


(n+t)At 

+ J exp [A ((n + 1) At - 1)] Bu (x) dt (17) 

nAt 


where At is the step size. The simplest approximation for the 
integral is to assume that u(x) is constant, u(t) - u(nAt). A 
better approximation is to assume that u varies linearly from u n 
to u"+i, estimating u n+1 as un + (u n - u n-1 ). The resulting 
algorithm is 


n+1 

x t 


= d> x" + 0 B (3u n - u"* 1 ) H 


(18) 


where e is the integral of the state-transition matrix 
Reference (11) describes a comparative evaluation of seven 
alternative structural integration algorithms including that of 
Eq. ( 18 ). The modified state-transition matrix integrator of Eq. 
(18) was shown to be superior to the others in terms of 
numerical stability and accuracy. 

For aeroelastic analysis, two steps are generally required in 
performing the calculations. In the first step, the steady-state 
flow field is calculated to account for wing thickness, camber, 
and mean angle of attack thus providing the starting flow field 
for the aeroelastic analysis. The second step is to prescribe an 
initial disturbance to begin the structural integration. 
Disturbance velocities in one or more modes, rather than 
displacements, have been found to be distinctly superior in 
avoiding nonphysical, strictly numerical transients and their 
possible associated instabilities. For the applications presented 
below, 1000 time steps were typically used to compute about 
three cycles of the dominant flutter mode and about 20 cycles of 
the higher frequency fourth mode (second torsion). In 
determining a flutter point, the freestream Mach number M and 
the associated freestream speed U were held fixed. A judiciously 
chosen value of the dynamic pressure pU2/2 is used to compute 
the free decay transients. These resulting transients of the 


generalized coordinates are analyzed (see below) for their 
content Of (tamped or growing sine-waves, with the rates of 
growth or decay indicating whether the dynamic pressure is 
above or below the flutter value. This analysis then indicates 
whether to increase or decrease the value of dynamic pressure 
in subsequent runs to determine a neutrally stable result. 
Further details are given in the following section on modal 
identification. 


Amplitude 



(a) aeroelastic transient and least-squares curve fit. 



(b) component modes from curve fit. 




« 


Fig. 1 Sample modal identification from free decay transient for 
a 45o sweptback wing calculated using CAP-TSD. 
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As previously mentioned, CAP-TSD generates free decay 
transients that must be analyzed for the modal stability 
characteristics. An example transient for a 45° sweptback 
wing, calculated using CAP-TSD is shown in Fig. 1(a). All four 
modes used in the analysis were excited by a velocity initial 
condition which produces a complicated decay record. This 
record is analyzed using a least-squares curve-fit with complex 
exponential functions with a program that is a derivative of the 
one described in Ref. (12). The components of the transient of 
Fig. 1(a) are plotted in Fig. 1(b) to the same scale. The free 
decay properties of each mode for this condition are readily 
apparent. A sufficient range of dynamic pressure must be 
considered to determine all relevant flutter points. 


WING FLUTTER TEST CASE 

To assess the CAP-TSD code for flutter applications, a simple 
well-defined wing case was selected as a first step toward 
performing aeroelastic analyses for complete aircraft 
configurations. The wing being analyzed is a semispan wind- 
tunnel-wall-mounted model that has a quarter-chord sweep 
angle of 45°, a panel aspect ratio of 1.65, and a taper ratio of 
0.66 (8). The wing is a proposed AGARD standard aeroelastic 
configuration (13) which was tested in the Transonic Dynamics 
Tunnel (TDT) at NASA Langley Research Center. A planview of 
the wing is shown in Fig. 2. The wing has a NACA 65A004 
airfoil section and was constructed of laminated mahogany. In 
order to obtain flutter for a wide range of Mach number and 
density conditions in the TDT, holes were drilled through the 
wing to reduce its stiffness. To maintain the aerodynamic shape 
of the wing, the holes were filled with a rigid foam plastic. A 
photograph of the model mounted in the TDT is shown in Fig. 3. 
The wing is being modeled structurally using the first four 
natural vibration modes which are illustrated in Figs. 4 and 5. 
Figure 4 shows oblique projections of the natural modes while 
Fig. 5 shows the corresponding deflection contours. These modes 
which are numbered 1 through 4 represent first bending, first 
torsion, second bending, and second torsion, respectively, as 
determined by a finite element analysis. The modes have natural 
frequencies which range from 9.6 Hz for the first bending mode 
to 91.54 Hz for the second torsion mode. 



Fig. 2 Planview of 45° sweptback wing. 


RESULTS AND DISCUSSION 

Results are presented from CAP-TSD to assess the code for 
flutter prediction. These initial results are evaluated by making 
comparisons with linear theory calculations and with the 
experimental flutter data of Ref. (8). 


Generalized Aerodynamic Forces (GAFs) 

Generalized aerodynamic forces were obtained using CAP- 
TSD to assess the accuracy of the code for applications with 
modal oscillations. Parallel calculations were performed using 
the FAST (14) linear-theory subsonic kernel-function program 
for comparison. For consistency of comparison, the linear 
potential equation (F * G - H « 0) option within CAP-TSD was 
selected and the wing was modeled as a flat plate (zero 
thickness). The results from CAP-TSD were obtained using the 
pulse transfer-function analysis (15). In the pulse analysis, 
the GAFs are computed indirectly from the response of the flow 
field due to a smoothly varying exponentially shaped pulse. A 
small pulse is prescribed in a given vibration mode and the 
aerodynamic transients are computed. The GAFs in the frequency 
domain are then determined by a transfer-function analysis 
involving Fast Fourier Transforms. This capability of CAP-TSD 
was recently developed and applied to study trends in 
aerodynamic forces (15). With the pulse analysis, the GAFs for 
one column of the aerodynamic matrix can be generated in a 
single run. A complete set of GAFs requires as many computer 
runs as the number of modes. 

Generalized aerodynamic force coefficients for the 45° 
sweptback wing are given in Fig. 6. These GAFs Ajj are defined 
as the force coefficients from the pressure induced by mode j 
acting through the displacements of mode i. The values 
represent work divided by the dynamic pressure and by c 2 . The 
results are subsequently presented in the form of real and 
imaginary components of An, A21, A12, and A22 as functions of 
reduced frequency k. The results from harmonic calculations 
from the FAST program are also compared with the CAP-TSD 
results. Good agreement is shown particularly in the A11 and 
A22 coefficients which are the first two diagonal terms of the 
aerodynamic matrix. The largest differences between the two 
sets of results occur for k > 0.5. The accuracy of the GAFs is 



Fig. 3 45° sweptback wing in the NASA Langley Transonic 

Dynamics Tunnel. 
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most important, however, for 0.1 < k < 0.5, since this is the 
range of reduced frequency where flutter typically occurs. The 
higher two modes, second bending and second torsion, have only a 
small influence on flutter for this case. Although there are 
slightly larger differences between results for the higher modes 
(not shown), the GAFs from CAP-TSD agree well with those 
from FAST. This good agreement thus verifies the CAP-TSD code 
for generalized force computation and indicates that the finite- 
difference grid is adequate for such applications. 

F lut te r J temlte 

Flutter calculations were performed for the 45° sweptback 
wing using CAP-TSD to assess the code for aeroelastic 
applications. Two sets of results are presented corresponding to: 


(1) using the linear potential equation and modeling the wing 
aerodynamlcally as a flat plate (zero thickness) and (2) using 
the complete (nonlinear) TSD equation and including wing 
thickness. The first set of results allows for direct comparison 
with parallel linear theory calculations performed using FAST. 
The second set of results more accurately models the wing 
geometry as well as the flow physics. All of the results are 
compared with the experimental flutter data of Ref. (8) which 
spans the range 0.338 * M 1 .141 . 

Comparisons of flutter characteristics from the linear 
calculations with the experimental data are given in Fig. 7. 
Plots of flutter speed index (defined as U/(bp o a Vp)) and 

nondimensional flutter frequency (defined as <n/a>a) as functions 






(b) Mode 2, f 2 -38.17Hz 



(d) Mode 4, f 4 -91.54 Hz 


Fig. 4 Oblique projections of natural vibration modes of 45° 
sweptback wing. 
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Fig. 5 Deflection contours of natural vibration modes of 45° 
sweptback wing. 


of freestream Mach number, are shown in Figs. 7(a) and 7(b), 
respectively. The experimental flutter data defines a typical 
transonic flutter "dip" with the bottom near M - 1.0 for this 
case. (Note that these results are shown with an expanded 
scale.) The bottom of the dip in flutter speed index (Fig. 7(a)) 
was defined by the approach to the M - 1.072 flutter point 
during the wind tunnel operation. Results from the CAP-TSD 
(linear) code are presented at twelve values of M covering the 
entire Mach number range over which the flutter data was 
measured. Results from the FAST program are presented for the 
limited range 0.338 0.96 since the method Is restricted to 

subsonic freestreams. Overall, the linear CAP-TSO results 
compare well with the experimental data for subsonic as well as 
supersonic Mach numbers. Note that the subsonic FAST results 
are also in good agreement with the data. Such a result is not 



(a) GAF resulting from pressure Induced by first 
bending acting through the first bending 
displacements. 



( b ) GAF resulting from pressures induced by first 
bending acting through the first torsion 
displacements. 


unexpected for this very thin wing of moderate sweep and taper 
at zero angle of attack. It does indicate that the wing properties 
are well-defined for benchmark purposes. 

In the subsonic Mach number range, the CAP-TSD and FAST 
calculations predict a slightly unconservative flutter speed, 
except at M - 0.338, by as much as 2% (Fig. 7(a)), and a 
higher flutter frequency (Fig. 7(b)) in comparison with the 
experimental data. In general though, the linear CAP-TSD 
results agree well with the FAST results in both flutter speed 
and frequency. The good agreement in this three-way 
correlation between experiment, linear theory, and CFD flutter 
results gives confidence in the CAP-TSD code for flutter 
prediction. 


■ CAP-TSD (linear) 



(c) GAF resulting from pressure induced by first 
torsion acting through the first bending 
displacements. 



(d) GAF resulting from pressure induced by first 
torsion acting through the first torsion 
displacements. 


Fig. 6 Comparisons between generalized aerodynamic forces 
(GAFs) calculated by CAP-TSD and FAST for the 45° 
sweptback wing at M ■ 0.9. 
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Comparisons of flutter characteristics from the linear and 
nonlinear CAP-TSD calculations with the experimental data are 
given in Fig. 6. Figure 8(a) shows flutter speed index versus 
Mach number and Fig. 8(b) shows nondimenslonal flutter 
frequency versus Mach number. Three flutter points are plotted 
from the nonlinear CAP-TSD calculations corresponding to 
M - 0.678, 0.901 , and 0.96. Comparisons between the two sets 
of CAP-TSD results show differences due to wing thickness and 
nonlinear effects. With increasing Mach number these 
differences become larger. For example, at M * 0.678, 0.901, 
and 0.96, the flutter speed index decreased by 1%, 5%, and 
19%, respectively, as shown in Fig. 8(a). Similar decreases 



0 .2 .4 .6 .8 1.0 1.2 1.4 

M 


( a ) flutter speed index versus Mach number. 



0 .2 .4 .6 .8 1.0 1.2 1.4 


( b ) nondimenslonal flutter frequency versus Mach 
number. 

Fig. 7 Comparisons between linear flutter calculations with 
experimental data for the 45° sweptback wing. 


also occur in the flutter frequency (Fig. 8(b)). The decrease in 
flutter speed at M - 0.901 is largely due to Including wing 
thickness since there are no supersonic points in the flow at this 
condition. The decrease in flutter speed at M - 0.96 is 
attributed to both wing thickness and nonlinear effects since an 
embedded supersonic region of moderate size was detected in the 
wing tip region. The nonlinear CAP-TSD results at both 
M * 0.901 and 0.96 are slightly conservative in comparison 
with the experimental flutter speed index value. Nonetheless, 
the nonlinear CAP-TSD flutter results compare favorably with 
the experimental data, which is the first step toward validating 
the code for general transonic aeroelastic applications. 



0 .2 .4 .6 .8 1.0 1.2 1.4 

M 


(a) flutter speed index versus Mach number. 



0 .2 .4 .6 .8 1.0 1.2 1.4 

M 


(b) nondimenslonal flutter frequency versus Mach 
number. 

Fig. 8 Comparisons between linear and nonlinear CAP-TSD 
flutter predictions with experimental data for the 45° 
sweptback wing. 
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CONCLUDING REMARKS 

The application and assessment ot the recently-developed 
CAP-TSD (Computational Aeroelasticity £rogram - Xransonic 
Small Disturbance) code for flutter prediction was presented 
and discussed. The CAP-TSD code has been developed for 
aeroelastic analysis of complete aircraft configurations in the 
flutter-critical transonic speed range. The code was previously 
applied to the calculation of steady and unsteady pressures on 
• wings and configurations at subsonic, transonic, and supersonic 
Mach numbers. Comparisons of these results with other 
methods and with experimental data have been favorable. 
A However, the purpose of CAP-TSD is for aeroelastic analysis. 
The evaluation and assessment of the CAP-TSD aeroelastic 
capability was the subject of the present study. Although the 
emphasis of the discussion was on CAP-TSD, the general 
computational procedures regarding time-marching flutter 
analysis apply to other computational fluid dynamics codes as 
well. 

Generalized aerodynamic forces (GAFs) and flutter 
boundaries were presented for a 45° sweptback wing. The GAFs 
from CAP-TSD, calculated using the linear potential equation, 
agreed well with those from a linear-theory subsonic kernel- 
function program. This agreement thus verifies CAP-TSD for 
generalized force computation. The flutter boundaries from 
CAP-TSD (linear) were in agreement with parallel subsonic 
linear theory results and compared well with the experimental 
flutter data for subsonic and supersonic freestream Mach 
numbers. The nonlinear CAP-TSD flutter results also compared 
favorably with the experimental data which is the first step 
toward validating the code for general transonic aeroelastic 
applications. 
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